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We present the design, fabrication, and characterization of a planar silicon photonic crystal 
cavity in which large position-squared optomechanical coupling is realized. The device consists 
of a double-slotted photonic crystal structure in which motion of a central beam mode couples 
to two high-Q optical modes localized around each slot. Electrostatic tuning of the structure is 
used to controllably hybridize the optical modes into supermodes which couple in a quadratic 
fashion to the motion of the beam. From independent measurements of the anti-crossing of 
the optical modes and of the optical spring effect, the position-squared vacuum coupling rate is 
measured to be as large as g /2n — 245 Hz to the fundamental in-plane mechanical resonance of the 
structure at cj m / 2tt = 8.7 MHz, which in displacement units corresponds to a coupling coefficient 
of g' / 27 t = 1 THz/nm 2 . This level of position-squared coupling is approximately five orders of 
magnitude larger than in conventional Fabry-Perot cavity systems. 


PACS numbers: 42.50.Wk, 42.65.-k, 62.25.-g 

I. INTRODUCTION 

In a cavity-optomechanical system the electromagnetic 
field of a resonant optical cavity or electrical circuit is 
coupled to the macroscopic motional degrees of freedom 
of a mechanical structure through radiation pressure pQ. 
Cavity-optomechanical systems come in a multitude 
of different sizes and geometries, from cold atomic 
gases [2| and nanoscale photonic structures |3|, to the 
kilogram/kilometer scale interferometers developed for 
gravitational wave detection [4]. Recent technological 
advancements in the field have led to the demonstration 
of optomechanically induced transparency 0 i, 
back-action cooling of a mechanical mode to its quantum 
ground state 00 , and ponderomotive squeezing of the 
light field [TOj [IT]. 

The interaction between light and mechanics in a 
cavity-optomechanical system is termed dispersive when 
it couples the frequency of the cavity to the position or 
amplitude of mechanical motion. To lowest order this 
coupling is linear in mechanical displacement, however, 
the overall radiation pressure interaction is inherently 
nonlinear due to the dependence on optical intensity. To 
date, this nonlinear interaction has been too weak to 
observe at the quantum level in all but the ultra-light 
cold atomic gases [2], and typically a large optical drive 
is used to parametrically enhance the optomechanical 
interaction. Qualitatively novel quantum effects are 
expected when one takes a step beyond the standard 
linear coupling and exploits higher order dispersive 
optomechanical coupling. In particular, “x 2 -coupling” 
where the cavity frequency is coupled to the square 
of the mechanical displacement has been proposed as 
a means for realizing quantum non-demolition (QND) 
measurements of phonon number HMD , measurement 


of phonon shot noise m, and the cooling and squeezing 
of mechanical motion p~6HT8] . In addition to dispersive 
coupling, an effective x 2 -coupling via optical homodyne 
measurement has also been proposed, with the capability 
of generating and detecting non-Gaussian motional 

states m- 

The dispersive x 2 -coupling between optical and 
mechanical resonator modes in a cavity-optomechanical 
system is described by the coefficient g' = d 2 uj c /dx 2 , 
where uj c is the frequency of the optical resonance of 
interest and x is the generalized amplitude coordinate 
of the displacement field of the mechanical resonance. 
One can show via second-order perturbation theory [2T3 
l2Tj that x 2 -coupling arises due to linear cross-coupling 
between the optical mode of interest and other modes 
of the cavity. In the case of two nearby resonant 
modes, the magnitude of the x 2 -coupling coefficient 
depends upon the square of the magnitude of the linear 
cross-coupling between the two modes ( g ) and inversely 
on their frequency separation or tunnel coupling rate 
(2J), g' = g 2 /2J. In pioneering work by Thompson, 
et al. 12], a Fabry-Perot cavity with an optically-thin 
SisN 4 membrane positioned in between the two end 
mirrors was used to realize x 2 -coupling via hybridization 
of the degenerate modes of optical cavities formed on 
either side of the partially reflecting membrane. More 
recently, a number of cavity-optomechanical systems 
displaying x 2 -coupling have been explored, including 
double microdisk resonators [22] . micro disk-cantilever 
systems m, microsphere-nanostring systems EU, 
atomic gases trapped in Fabry-Perot cavities [2], and 
paddle nano-cavities [211 . 

Despite significant technical advances made in recent 
years OEaEaES], the use of x 2 -coupling for measuring 
or preparing non-classical quantum states of a mesoscopic 
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mechanical resonator remains an elusive goal. This is a 
direct result of the small coupling rate to motion at the 
quantum level, which for x 2 -coupling scales as the square 
of the zero-point motion amplitude of the mechanical 
resonator, x 2 pf = where m is the motional mass 

of the resonator and c is the resonant frequency. As 
described in Ref. [14], one method to greatly enhance 
the x 2 -coupling in a multi-mode cavity-optomechanical 
system is to fine tune the mode splitting 2J to that of 
the mechanical resonance frequency. 

In this work we utilize a quasi two-dimensional 
(2D) photonic crystal structure to create an optical 
cavity supporting a pair of high-Q optical resonances 
in the 1500 nm wavelength band exhibiting large linear 
optomechanical coupling. The double-slotted structure 
is split into two outer slabs and a central nanobeam, 
all three of which are free to move, and electrostatic 
actuators are integrated into the outer slabs to allow for 
both the trimming of the optical modes into resonance 
and tuning of the tunnel coupling rate J. Due to 
the form of the underlying photonic bandstructure the 
spectral ordering of the cavity supermodes in this 
structure may be reversed, enabling arbitrarily small 
values of J to be realized. Measurement of the optical 
resonance anti-crossing curve, along with calibration of 
the linear optomechanical coupling through measurement 
of the dynamic optical spring effect, yields an estimated 
x 2 -coupling coefficient as large as g r /2n = 1 THz/nm 2 
(g'/2 tt = 245 Hz) to the fundamental mechanical 
resonance of the central beam at cj m /27r = 8.7 MHz. 
Additional measurements of g' through the dynamic 
and static optical spring effects are also presented. 
This level of x 2 -coupling is approximately five orders 
of magnitude larger than in conventional Fabry-Perot 
MIM systems demonstrated to date [26], and two orders 
of magnitude larger than in the smaller mode volume 
fiber-gap cavities [25] . 


II. THEORETICAL BACKGROUND 

Before we discuss the specific double-slotted photonic 
crystal cavity-optomechanical system studied in this 
work, we consider a more generic multi-moded system 
consisting of two optical modes which are dispersively 
coupled to the same mechanical mode, and in which the 
dispersion of each mode is linear with the amplitude 
coordinate x of the mechanical mode. If we further 
assume a purely optical coupling between the two 
optical modes, the Hamiltonian for such a three-mode 
optomechanical system in the absence of drive and 
dissipation is given by H = Hq + Hqm + Hj\ 


Hq = fyjj\a\a\ + f\juj2Ci\a2 + ( 1 ) 

#om = h(gia\ai + g 2 a\a 2 )x, (2) 

H j = hJ (a\a 2 + a\ai). (3) 



Figure 1. (a) Double-slotted photonic crystal cavity with 

optical cavity resonances ( 01 , 02 ) centered around the two 
slots, and three fundamental in-plane mechanical resonances 
corresponding to motion of the outer slabs ( 61 , 62 ) and the 
central nanobeam (63). Tuning the equilibrium position of the 
outer slabs 61 and 62 , and consequently the slot size on either 
side of the central nanobeam, is achieved by pulling on the 
slabs (red arrows) through an electrostatic force proportional 
to the square of the voltage applied to capacitors on the outer 
edge of each slab, (b) Dispersion of the optical modes as 
a function of X 3 , the in-plane displacement of the central 
nanobeam from its symmetric equilibrium position. Due to 
tunnel coupling at a rate J the slot modes a\ and <22 hybridize 
into the even and odd supermodes a+ and a_, which have 
a parabolic dispersion near the central anti-crossing point 
(ui as CJ 2 ). (c) SEM image of a fabricated double-slotted 

photonic crystal device in the SOI material system, (d) 
Zoom-in SEM image showing the capacitor gap (~ 100 nm) 
for the capacitor of one of the outer slabs, (e) Zoom-in SEM 
image showing some of the suspending tethers of the outer 
slabs which are of length 2.5 gm and width 155 nm. The 
central beam, which is much wider, is also shown in this 
image. 


Here, and Ui are the annihilation operator and the 
bare resonance frequency of the ith optical resonance, 
x = (6^+6)x zp f is the quantized amplitude of motion, x zp f 
the zero point amplitude of the mechanical resonance, 
is the bare mechanical resonance frequency, and gi is the 
linear optomechanical coupling constant of the ith optical 
mode to the mechanical resonance. Without loss of 
generality, we take the bare optical resonance frequencies 
to be equal (ui = 0 J 2 m cjo), allowing us to rewrite the 
Hamiltonian in the normal mode basis a± = (di ± 62 )/ y/2 
as, 
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H = 7kj+(0)a+a+ + /kj_(0)ala_ + hujm&b 

+ h ^ ^ + ala_^ x 

+ h ^ (®+“- + «ia + ^ 5, (4) 

where u;±(0) = cjo ± J. 

For |J| >> u; m such that x can be treated as a 
quasi-static variable na in, the Hamiltonian can be 
diagonalized resulting in eigenfrequencies u±(x), 


/->, I (dl + 92 ) ~ | J f(dl — 92) 2 /r\ 

CJ±(x)^CJ 0 + - 2 - xd = J ( 1 + -^J2- x )• ( 5 ) 

As shown below, in the case of the fundamental in-plane 
motion of the outer slabs of the double-slotted photonic 
crystal cavity we have only one of g\ or c /2 non-zero, 
whereas in the case of the fundamental in-plane motion 
of the central nanobeam we have g\ ~ — # 2 - 

For a system in which the mechanical mode couples 
to the a\ and a 2 optical modes with linear dispersive 
coupling of equal magnitude but opposite sign (gi = 
—#2 = g), the dispersion in the quasi-static normal 
mode basis is purely quadratic with effective x 2 -coupling 
coefficient, 


9' = 9 2 /2J, (6) 

and quasi-static Hamiltonian, 

H « h ((j + (0) + g'x 2 ) n + 

+ h (u~ (0) - g'x 2 ) n_ + fvu m n b , (7) 

where h± are the number operators for the a± 
supermodes and h & is the number operator for the 
mechanical mode. Rearranging this equation slightly 
highlights the interpretation of the x 2 optomechanical 
coupling as inducing a static optical spring, 


where g' = g'x 2 pf = g 2 / 2J and g = gx zp f are the x 2 
and linear vacuum coupling rates, respectively. It is 
tempting to assume from eq. § that by monitoring the 
optical transmission through the even or odd supermode 
resonances, that one can then perform a continuous 
quantum non-demolition (QND) measurement of the 
phonon number in the mechanical resonator QJE7H29]. 
As noted in Refs. esi m, however, the quasi-static 
picture described by the dispersion of eq. (J5| fails to 
capture residual effects resulting from the non-resonant 
scattering between the a + and a_ supermodes which 
depends linearly on x (last term of eq. Only in the 
vacuum strong coupling limit (g/ft > 1) can one realize 
a QND measurement of phonon number mm- 

The regime of \2J\ ~ uo m is also very interesting, 
and explored in depth in Refs. [I4j [30] . Transforming 
to a reference frame which removes in eq. 0 the 
radiation pressure interaction between the even and 
odd supermodes to first order in g , yields an effective 
Hamiltonian given by [I4j [3T, 


i7 e ff ~ /kj + (0)n_|_ + hoj_ (0)rz_ 

1 1 




, h g 2 


2 J ^rri 

1 

2 J — uj n 


2J T UJrH 

1 


a+a+ - dla_ 


2 J T uj r 


)(*+») 

^a+a_+aLa + j , (10) 


where we assume \g/S\ 1 for 4 = \2J\ — cj m , and terms 

of order g 3 /(2J ± c o m ) 2 and higher have been neglected. 
In the limit \J\^> we recover the quasi-static result 
of eq. 0 - whereas in the near-resonant limit of \5\ <C 
|J|,o; m we arrive at, 


H e ff « /kj + (0)n + + /kj_(0)n_ + hhJm^b 

Q 2 

+ [ 2 sgn(J) (n+ - n_) (n b + 1) 

+ 2 n+n_ + n++ n_]. (11) 


H ~ /kj + (0)fi + + ficj_(0)n_ 

+ ft [w m n 6 + g' (n+ -n_)£ 2 ] , (8) 

where the static optical spring constant k s = 
2hg'(n+— n_) depends upon the average intra-cavity 
photon number in the even and odd optical supermodes, 
n± = (n±). 

For a sideband resolved system (cj m /c), the 

quasi-static Hamiltonian can be further approximated 
using a rotating-wave approximation as, 

H ~ /i[cj+(0) + 2 g'(fib + l/2)]n_|_ 

+ fi[cj_(0) — 2 g' {fib + l/2)]n_ + (9) 


Here we have neglected highly oscillatory terms such as 
(a^d_) 2 and 6 2 , a good approximation in the sideband 
resolved regime (k, | J|). From eq. © we 

find that the frequency shift per phonon of the optical 
resonances is much larger than in the quasi-static case 
{g 2 /2\5\ g 2 /2\J\). Although a QND measurement of 
phonon number still requires the vacuum strong coupling 
limit, this enhanced read-out sensitivity is attainable 
even for gj k <C 1. Equation © also indicates that, 
much like the QND measurement of phonon number, in 
the near-resonant limit a measurement of the intra-cavity 
photon number stored in one optical supermode can 
be performed by monitoring the transmission of light 
through the other supermode EHIHI]- 
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III. DOUBLE-SLOTTED PHOTONIC CRYSTAL 
OPTOMECHANICAL CAVITY 

A sketch of the double-slotted photonic crystal cavity 
structure is shown in Fig. EK As detailed below 
and elsewhere [32], the optical cavity structure can be 
thought of as formed from two coupled photonic crystal 
waveguides, one around each of the nanoscale slots, and 
each with propagation direction along the z-axis. A 
small adjustment 5%) in the lattice constant is used 
to produce a local shift in the waveguide band-edge 
frequency, resulting in trapping of optical resonance to 
this “defect” region. Optical tunneling across the central 
photonic crystal beam, which in this case contains only 
a single row of holes, couples the cavity mode of slot 1 
(ai) to the cavity mode of slot 2 (( 22 ). 

The two outer photonic crystal slabs and the central 
nanobeam are all mechanically compliant, behaving as 
independent mechanical resonators. The mechanical 
resonances of interest in this work are the fundamental 
in-plane flexural modes of the top slab, the bottom slab, 
and the central nanobeam, denoted by 61 , 62 and 63 , 
respectively. For a perfectly symmetric structure about 
the z-axis of the central nanobeam, the linear dispersive 
coupling coefficients of the 63 mode of the central 
nanobeam to the two slot modes a\ and a 2 are equal in 
magnitude but opposite in sign, resulting in a vanishing 
linear coupling at the resonant point where uq = Coq (c.f., 
eq. <§)• Figure [TJd shows a plot of the dispersion of 
the optical resonances as a function of the nanobeam’s 
in-plane displacement (£ 3 ), illustrating how the linear 
dispersion of the slot modes (ai,a 2 ) transforms into 
quadratic dispersion for the upper and lower supermode 
branches (a + ,a_) in the presence of tunnel coupling J. 
The mechanical modes of the outer slabs (fq, 62 ) provide 
degrees of freedom for post-fabrication tuning of the 
slotted waveguide optical modes, i.e., to symmetrize the 
structure such that uj\ = uq. This is achieved in practice 
by integrating metallic electrodes which form capacitors 
at the outer edge of the two slabs of the structure as 
schematically shown in Fig. EK 

The double-slotted photonic crystal cavity of this 
work is realized in the silicon-on-insulator (SOI) material 
system, with a top silicon device layer thickness of 
220 nm and an underlying buried oxide (BOX) layer 
of 3 fi m. Fabrication begins with the patterning of the 
metal electrodes of the capacitors, and involves electron 
beam (ebeam) lithography followed by evaporation and 
lift-off of a bi-layer consisting of a 5 nm sticking layer of 
chromium and a 150 nm layer of gold. After lift-off we 
deposit uniformly a ~ 4 nm protective layer of silicon 
dioxide. A second electron beam lithography step is 
performed, aligned to the first, to form the pattern of the 
photonic crystal and the nanoscale slots which separate 
the central nanobeam from the outer slabs. At this 
step, we also pattern the support tethers of the outer 
slabs and the cut lines which define and isolate the outer 
capacitors. A fluorine based (C 4 F 8 and SFe) inductively 


coupled reactive-ion etch (ICP-RIE) is used to transfer 
the ebeam lithography pattern through the silicon device 
layer. The remaining ebeam resist is stripped using 
trichloroethylene, and then the sample is cleaned in a 
heated Piranha (H 2 S 04 :H 2 02 ) solution. The devices 
are then released using a hydrofluoric (HF)acid etch to 
remove the sacrificial BOX layer (this also removes the 
deposited protective silicon dioxide layer), followed by a 
water rinse and critical point drying. 

A scanning electron microscope (SEM) image showing 
the overall fabricated device structure is shown in Fig.JTJ^. 
Zoom-ins of the capacitor region of one of the outer 
slabs and the tether region at the end of the nanobeam 
are shown in Figs. Et and e, respectively. Note that 
the geometry of the capacitors and the stiffness of the 
support tethers determine how tunable the structure is 
under application of voltages to the capacitor electrodes. 
The outermost electrode of each slab is connected to 
an independent low-noise DC voltage source, while the 
innermost electrodes are connected to a common ground, 
thereby allowing one to independently pull on each outer 
slab with voltages Vi and V 2 . In this configuration, we 
are limited to increasing the slots defining the optical 
modes around the central nanobeam. 


A. Photonic bandstructure 

To further understand the optical properties of the 
double-slotted photonic crystal cavity, we display in 
Fig. §1 the photonic bandstructure of the periodic 
waveguide structure. The parameters of the waveguide 
are given in the caption of Fig. UK Here we only 
show photonic bands that are composed of waveguide 
modes with even vector symmetry around the “vertical” 
mirror plane (cr z ), where the vertical mirror plane is 
defined by the z-axis normal and lies in the middle of 
the thin-film silicon slab. The fundamental (lowest lying) 
optical waveguide bands are of predominantly transverse 
(in-plane) electric field polarization, and are thus called 
TE-like. In the case of a perfectly symmetric structure, 
we can further classify the waveguide bands by their 
odd or even symmetry about the “horizontal” mirror 
plane (<j y ) defined by the y -axis normal and cutting 
through the middle of the central nanobeam. The two 
waveguide bands of interest that lie within the quasi-2D 
photonic bandgap of the outer photonic crystal slabs, 
shown as bold red and black curves, are labeled “even” 
and “odd” depending on the spatial symmetry with 
respect to a y of their mode shape for the dominant 
electric field polarization in the ^/-direction, E y (note that 
this labeling is opposite to their vector symmetry). The 
E y spatial mode profiles at the X-point for the odd and 
even waveguide supermodes are shown in Figs. [ 2)3 and c, 
respectively. 

An optical cavity is defined by decreasing the lattice 
constant 4.5% below the nominal value of <20 = 480 nm 
for the middle five periods of the waveguide (see Fig. IP- 
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Figure 2. (a) Bandstructure diagram of the periodic (along 
x) double-slotted photonic crystal waveguide structure. Here 
we only show photonic bands that are composed of modes 
with even vector symmetry around the “vertical” ( cr z ) mirror 
plane. The two waveguide bands of interest lie within the 
quasi-2D photonic bandgap of the outer photonic crystal 
slabs and are shown as bold red and black curves. These 
waveguide bands are labeled “even” (bold black curve) and 
“odd” (bold red curve) due to the spatial symmetry of their 
mode shape for the dominant electric field polarization in 
the y-direction, E y . The simulated structure is defined by 
the lattice constant between nearest neighbor holes in the 
hexagonal lattice (ao = 480 nm), the thickness of the silicon 
slab (d — 220 nm), the width of the two slots (s = 100 nm), 
and the refractive index of the silicon layer (nsi = 3.42). 
The hole radius in the outer slabs and the central nanobeam 
is r — 144 nm. The grey shaded region represents a 
continuum of radiation modes which lie above the light cone 
for the air cladding which surrounds the undercut silicon slab 
structure, (b) Normalized E y field profile at the X-point of 
the odd waveguide supermode, shown for several unit cells 
along the x guiding axis, (c) E y field profile of the even 
waveguide supermode. Waveguide simulations of (a-c) were 
performed using the plane-wave mode solver MPB [331 El- 
Normalized E y field profile of the corresponding localized 
cavity supermodes of (d) odd and (e) even spatial symmetry 
about the horizontal mirror plane. The lattice constant ao 
is decreased by 4.5% for the central five lattice constants 
between the dashed lines to localize the waveguide modes. 
Simulations of the full cavity modes were performed using the 
COMSOL finite-element method mode solver package m 


This has the effect of locally pushing the bands toward 
higher frequencies [Ml 132, which creates an effective 
potential that localizes the optical waveguide modes 
along the x-axis of the waveguide. The resulting 
odd and even TE-like cavity supermodes are shown 
in Figs. and e, respectively. These optical modes 
correspond to the normal modes a + and a_ in Section [TH 
which are symmetric and anti-symmetric superpositions, 
respectively, of the cavity modes localized around each 
slot (ai and < 12 ). Due to the non-monotonic decrease in 
the even waveguide supermode as one moves away from 
the X-bandedge (c.f., Fig. [2 ^l) , we find that the simulated 
optical Q-factor of the even a + cavity supermode is 
significantly lower than that of the odd a_ cavity 


Figure 3. Symmetric tuning of the slot widths of the 
double-slotted photonic crystal cavity showing (a) the mean 
wavelength shift and (b) the splitting 2J = — uj- of the 

even and odd cavity supermodes versus slot width s = s\ — 
S 2 . (c-d) Avoided crossing of the cavity supermodes obtained 
by tuning s 1 while keeping S 2 fixed at (c) S 2 = 90 nm and (d) 
S 2 = 95 nm. For all simulations in (a-d) the parameters of the 
cavity structure are the same as in Fig. [2] except for the slot 
widths. The simulations were performed using the COMSOL 
FEM mode solver [35] . 


supermode. This will be a key distinguishing feature 
found in the measured devices as well. 


B. Optical tuning simulations 

The slot width in the simulated waveguide and cavity 
structures of Fig. [2] is set at s = 100 nm. For this 
slot width we find a lower frequency for the even (a+) 
supermode than the odd (a_) supermode at the X-point 
photonic bandedge of the periodic waveguide and in the 
case of the localized cavity modes. Figure [3] presents 
finite-element method (FEM) simulations of the optical 
cavity for slot sizes swept from 90 nm to 100 nm in steps 
of 1 nm, all other parameters the same as in Fig. [2j For 
the slot widths tuned symmetrically (si = $2 = s), the 
mean wavelength of the even and odd cavity supermodes 
and their frequency splitting 2 J = — uj- are plotted 

in Fig. and Fig. respectively. As expected 
the mean wavelength drops for increasing slot width. 
The frequency splitting, however, also monotonically 
decreases with slot width, going from a positive value 
for s = 90 nm to a negative for s = 100 nm slots and 
crossing zero for a slot width of s = 95 nm. In Figs. IH 
and d the symmetry is broken by keeping 52 fixed and 
scanning .si; the cavity supermodes are driven through 
an anti-crossing with a splitting determined by the fixed 
slot width 52 . 

The spectral inversion of the even a + and odd 
a_ cavity supermodes predicted in Fig. originates 
in the unequal overlap of each mode with the air 
slots separating the two outer slabs from the central 
nanobeam. The odd supermode tends to be pushed 
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further from the middle of the central nanobeam, having 
slightly larger overlap with the air slots. An increase 
in the air region for increased slot size leads to a blue 
shift of both cavity supermodes. The odd mode having 
a larger electric field energy density in the air slots 
than the even mode is more affected by a change in 
the slot widths. Therefore, upon equal increase of the 
slot widths, the odd mode experiences larger frequency 
shifts than the even mode, which results in a tuning 
of the frequency splitting. For particular geometrical 
parameters of the central nanobeam [32], a change in the 
slot widths is sufficient to invert the spectral ordering 
of the supermodes. This means that arbitrarily small 
splittings can potentially be realized, which is important 
for applications in x 2 detection where the splitting enters 
inversely in the coupling (for the quasi-static case). 


IV. EXPERIMENTAL MEASUREMENTS 

Optical testing of the fabricated devices is performed 
in a nitrogen-purged enclosure at room temperature and 
pressure. A dimpled optical fiber taper is used to locally 
excite and collect light from the photonic crystal cavity, 
details of which can be found in Ref. [38]. The light from 
a tunable, narrow-bandwidth laser source in the telecom 
1550 nm wavelength band (New Focus, Velocity series) is 
evanescently coupled from the fiber taper into the device 
with the fiber taper guiding axis parallel with that of 
the photonic crystal waveguide axis, and the fiber taper 
positioned laterally at the center of the nanobeam and 
vertically a few hundreds of nanometers above the surface 
of the silicon chip. Relative positioning of the fiber taper 
to the chip is accomplished using a multi-axis set of 
encoded DC-motor stages with 50 nm step resolution. 
The polarization of the light in the fiber is polarized 
parallel with the surface chip in order to optimize the 
coupling to the in-plane polarization of the cavity modes. 

With the taper placed suitably close to a photonic 
crystal cavity (^ 200 nm), the transmission spectrum 
of the laser probe through the device features resonance 
dips at the supermode resonance frequencies, as 
shown in the intensity plots of Figs. [4^i- c. The 
resonance frequencies of the cavity modes are tuned 
via displacement of the top and bottom photonic 
crystal slabs, which can be actuated independently using 
their respective capacitor voltages V\ and W The 
capacitive force is proportional to the applied voltage 
squared m, and thus increasing the voltage V on 
a given capacitor widens the waveguide slot and 
(predominantly) increases the slot mode frequency cq 
(note the other optical slot mode frequency also increases 
slightly). For the devices studied in this work, the slab 
tuning coefficient with applied voltage (<a cap ) is estimated 
from SEM analysis of the resulting structure dimensions 
and FEM electromechanical simulations to be <a cap = 
25 pm/V 2 . 

We fabricated devices with slot widths targeted 


for a range of 75-85 nm, chosen smaller than the 
expected zero-splitting slot width of s = 95 nm so 
that the capacitors could be used to tune through 
the zero-splitting point. While splittings larger than 
150 GHz were observed in the nominal 85 nm slot width 
devices, splittings as small as 10 GHz could be resolved 
in the smaller 75 nm slot devices. As such, in what 
follows we focus on the results from a single device with 
as-fabricated slot size of s « 75 nm. 


A. Anti-crossing measurements 

Figure [4] shows intensity plots of the normalized 
optical transmission through the optical fiber taper when 
evanescently coupled to the photonic crystal cavity of 
a device with nominal slot width s = 75 nm. Here a 
series of optical transmission spectrum are measured by 
sweeping the probe laser frequency and the voltage Vf, 
with V 2 fixed at three different values. The estimated 
anti-crossing splitting from the measured dispersion of 
the cavity supermodes is 2J/27T = 50 GHz, 12 GHz, and 
—25 GHz for V 2 = 1 V, 15 V, and 18 V, respectively. 
In order to distinguish between the odd and even cavity 
supermodes at the anti-crossing point, we use the fact 
that both the coupling rate to the fiber taper K e and 
the intrinsic linewidth depend upon the symmetry 
of the cavity mode. First, the odd supermode branch 
becomes dark at the anti-crossing because it cannot 
couple to the symmetric fiber taper mode. Second, in 
the vicinity of the anti-crossing point the linewidth of 
the odd supermode branch narrows while the linewidth 
of the even supermode branch broadens [32]. Far from 
the anti-crossing region, the branches are asymptotic to 
individual slot modes and their linewidths and couplings 
to the fiber taper are similar. 

These features are clearly evident in the optical 
transmission spectra of Figs. & -c, as well as in the 
measured linewidth of the optical supermode resonances 
shown in Figs. [4^-h. Figure [4^i was taken with a small 
voltage V 2 = 1 V, corresponding to a small slot width 
at the anti-crossing point, and thus consistent with the 
even mode frequency being higher than the odd mode 
frequency for small slot widths (c.f., Fig. [3 ]d) . The exact 
opposite identification is made in Fig. H where V 2 = 18 V 
is much larger, corresponding to a larger slot width at the 
anti-crossing point. Fig. 0 > with V 2 = 15 V is close to the 
zero-splitting condition. For comparison, a simulation of 
the expected anti-crossing curves are shown in Figs. [4]i, 
e, and f for 82 = 93, 95, and 97 nm, respectively. Here we 
have taken the even superposition of the slot modes to 
have a lower Q-factor than the odd superposition of the 
slot modes, and the coupling of the fiber taper to be much 
stronger to the even mode than the odd mode, consistent 
with results from numerical FEM simulations. Good 
qualitative correspondence is found with the measured 
transmission curves of Figs. |§t-c. 

An estimate of the x 2 -coupling coefficient g' b3 can found 
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Figure 4. (a-c) Optical transmission measurements versus the wavelength of the probe laser showing the cavity mode 

anti-crossing and tuning of the photon tunneling rate. In these measurements the probe laser wavelength (horizontal axis) 
is scanned across the optical cavity resonances as the voltage across the first capacitor Vi is swept from low to high (vertical 
axis shows V 2 in V 2 , proportional to slab displacement). The second capacitor is held fixed at (a) V 2 = IV, (b) V 2 = 15V and 
(c) V 2 = 18V. The colorscale indicates the fractional change in the optical transmission level, AT, with blue corresponding 
to AT = 0 and red corresponding to AT « 0.25. From the three anti-crossing curves we measure a splitting 2 J equal to (a) 
50 GHz, (b) 12 GHz, and (c) —25 GHz. (d-f) Corresponding simulations of the normalized optical transmission spectra for the 
slot width si varied and the second slot width held fixed at (d) S 2 = 93 nm, (e) s 2 — 95 nm and (f) S 2 = 97 nm. The dispersion 
and tunneling rate of the slot modes are taken from simulations similar to that found in Fig. [ 3 ] (g) and (h) show the measured 
linewidths of the high frequency upper (black) and low frequency lower (red) optical resonance branches as a function of V?, 
extracted from (a) and (c), respectively. The narrowing (broadening) is a characteristic of the odd (even) nature of the cavity 
supermode, indicating the inversion of the even and odd supermodes for the two voltage conditions V 2 = 1 V and V 2 = 18 V. 
The lines are guides for the eye. 


from the simulated value of a cap and a fit to the measured 
tuning curves of Fig.[4]away from the anti-crossing point. 
Consider the anti-crossing curve of Fig. with the 
smallest splitting. Far from the anti-crossing point the 
tuning of the a\ and a 2 slot modes are measured to be 
linear with the square of V\\ g ai y? /27T = 3.9 GHz/V 2 
and g a2 ,v j 2 /2tt = 0.5 GHz/V 2 . For the simulated value of 
c^cap = 0.025 nm/V 2 the corresponding linear dispersive 
coefficients versus the first slot width are g aijl 5 Sl /27r = 
156 GHz/nm and g a 2 : s Sl /^ = 20 GHz/nm. Noting 
that a displacement amplitude x% for the fundamental 
in-plane mechanical mode of the central nanobeam is 
approximately equivalent to a reduction in the width 
of one slot by -X 3 and an increase in the other slot 
by +£ 3 , the linear optomechanical coupling coefficient 
between optical slot mode a\ and mechanical mode 
b 3 is estimated to be g aiM « (g ai ,s Sl + 9 ai ,-Ss 2 ) = 
(g ai ,6si -Poz/si) = 27r[136 GHz/nm], where by symmetry 
ga 1 —Ss 2 = — 90 . 2 , 68 !- Along with a measured splitting of 
2J/27T = 12 GHz, this yields through eq. © an estimate 
for the x 2 -coupling coefficient of g* h J2n « 1.54 THz/nm 2 . 


B. Transduction of mechanical motion 


Figure. [5] shows the evolution of the 
optically-transduced mechanical noise power spectral 
density (PSD) near the anti-crossing region of Fig. [4^i. 
In this plot 52 is fixed and si is varied over an estimated 
range of Ss± = ±0.3 nm around the anti-crossing. Optical 
motion is imprinted as intensity modulations of the 
probe laser which is tuned to the blue side of the upper 
frequency supermode. Here we choose the detuning 
point corresponding to A^ = ujl — ~ «/2\/3, 

where ujl is the probe laser frequency and k is the 
full-width at half-maximum linewidth of the optical 
resonance. This detuning choice ensures (maximal) 
linear transduction of small fluctuations in the 
frequency of the cavity supermode, which allows us 
to relate nonlinear transduction of motion with true 
nonlinear optomechanical coupling mi m- a probe 
power of Pi n = 10 /xW is used in order to avoid 
any nonlinear effects due to optical absorption, and 
the transmitted light is first amplified through an 
erbium-doped fiber amplifier before being detected 
on a high gain photoreceiver (transimpedance gain 










mechanical frequency (MHz) 


Figure 5. RF photocurrent noise spectrum for the optically 
transmitted light past the double-slotted photonic crystal 
cavity. Here the applied voltage V 2 = 1 V is held fixed and Vi 
is swept from just below to just above the anti-crossing point 
of Fig. S*- In these measurements the probe laser power is 
10 fi W at the input to the cavity, the probe laser frequency 
is set on the blue side of the upper frequency supermode 
resonance, A l ~ «/2>/3, and the fiber taper is placed in 
the near-field of the photonic crystal cavity resulting in an 
on-resonance dip in transmission of approximately AT = 
15%. The vertical axis in this plot is converted to a change 
in the slot width, 5si, using the numerically simulated value 
of <a cap = 0.025 nm/V 2 . The color indicates the magnitude of 
the RF noise in dBm/Hz, where the colorscale from 0-14 MHz 
is shown on the left of the scalebar and the colorscale from 
14-20 MHz is shown on the right of the scalebar (a different 
scale is used to highlight the noise out at 2 uo b3 . 


of the bi mode with the lower frequency mechanical 
resonance is made possible due to the increasing signal 
transduction of this resonance as si is increased above 
the anti-crossing point. Since we are probing the upper 
frequency optical supermode, for s 1 > 52 (5s 1 > 0) 
the supermode is approximately a\ which is localized to 
slot 1 and sensitive primarily to the motion of b\. We 
see an opposite trend for the 62 resonance, with larger 
transduction gain for si < S 2 (5s 1 < 0). The frequencies 
of both these modes is lower than found in numerical 
simulations, likely due to squeeze-film damping effects 
not captured in the FEM analysis m ■ 

The noise peak at c<j m /27r = 8.73 MHz behaves 
altogether differently than the b\ and 62 resonances, and 
is identified with the 63 mode of the central nanobeam 
(although again at a lower frequency than expected from 
FEM simulation). This noise peak is transduced with 
roughly equal signal levels for 5s 1 > 0 and 5s 1 < 0 , 
but significantly drops in strength for ~ 0 near 
the anti-crossing. This is the expected characteristic of 
the 63 mode, where the dispersive linear optomechanical 
coupling to the 63 should vanish at the anti-crossing 
point. Also shown in Fig. [5^i is the noise at 2wb 3 /27r ~ 
17.5 MHz, which shows a weakly transduced resonance 
with signal strength peaked around 5s 1 = 0. The 
suppression in transduction of the noise peak at u ;& 3 
concurrent with the rise in transduction of the noise 
peak at 2 c ^ 3 is a direct manifestation of the transition 
from linear (g ai ,b 3 or g a2 ,b 3 ) to position-squared (g' bs ) 
optomechanical coupling. 


10 4 V/A, NEP= 12 pW/Hz 1 / 2 , bandwidth 150 MHz). 
The resulting radio-frequency (RF) photocurrent noise 
spectrum is plotted in Fig. [5| 

To help identify the measured noise peaks numerical 
FEM simulations of the mechanical properties of 
the double-slotted structure were performed. Taking 
structural dimensions from SEM images, the simulated 
mechanical frequency for the fundamental in-plane 
resonances of the two outer slabs (bi and 62 ) is found 
to be / 2tt = 8.4MHz. An effective motional mass 
for the slab modes of m = 35 pg was determined by 
integrating, over the volume of the structure, the mass 
density of the silicon slab weighted by the normalized, 
squared displacement amplitude of the slab’s motion [39] . 
The corresponding estimate of the zero-point amplitude 
of the slab modes is given by x zp f = (h/2muj rn ) 1 / 2 = 
5.6 fm. The resonance frequency, effective motional mass, 
and zero-point amplitude for the fundamental in-plane 
resonance of the central nanobeam ( 63 ) are simulated to 
be u; m /27r = 10.7MHz, m = 3.6 pg, and x zp f = 15.4 fm, 
respectively. 

Comparing to Fig. [5^i, the two lowest frequency noise 
peaks are thus identified as due to the thermal motion of 
the bi and 62 modes of the outer slabs, with uj bl /2ir = 
5.54 MHz and cj^ 2 /27r = 6.34 MHz. The identification 


C. Static and dynamic optical spring 
measurements 

Our previous estimate of g ' h3 from the anti-crossing 
curves relied on the approximate correspondence between 
the static displacement of the outer slabs and the 
fundamental in-plane vibrational amplitude of the 63 
mode of the central nanobeam. A more accurate 
determination of the true x 2 -coupling coefficient to 63 
can be determined from two different optical spring 
measurements. Far from the anti-crossing one can 
determine the linear optomechanical coupling coefficient 
between the optical slot modes and the 63 mechanical 
mode from the dynamic back-action of the intra-cavity 
light field on the mechanical frequency, which in 
conjunction with the measured anti-crossing splitting 
yields g ' b3 via eq. (J6|. A direct measurement of g bs can 
also be obtained from the static optical spring effect near 
the anti-crossing point as indicated in eq. (|9|. 

Figure [ 6^1 shows the dependence of the mechanical 
resonance frequency of the 63 mode of the central 
nanobeam versus the laser detuning A l when the device 
is tuned far from the anti-crossing point in Fig. [4^i 
(Vi = 1 V and V 2 = 1 V). In these measurements 
the probe laser power is fixed at P in = 10 fi W and 
the laser frequency is scanned across the upper optical 
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Figure 6 . (a) Dynamic optical spring effect measured by 

exciting the upper frequency supermode resonance far from 
the anti-crossing point (~ a± mode) [Vi = V 2 = 1 V, 
P in = 10 /xW, k/2it = 12.5 GHz, AT « 10%]. (b) Static 
optical spring shift of the 63 resonance frequency versus laser 
detuning A l from the upper (~ 0 +) supermode resonance 
near the anti-crossing point [Vi = 10.9 V, V 2 = IV, 
Pin = 50 /iW, k,/2tt = 26 GHz, AT ^ 25%]. In both (a) 
and (b) V 2 is fixed at 1 V (see Fig. Bo and the measured 
data (circles) correspond to a Lorentzian fit to the resonance 
freqency of the optically transduced thermal noise peak at 
ujb 3 . In (a) the red curve is a fit to the data using a dynamical 
optical spring model m with linear optomechanical coupling 
coefficient g ai ,b 3 / 27r s= 1.72 MHz. In (b) the red curve is a 
fit to the data using a static spring model (c.f., eqs. © and 
([ 9 ])) with £ 2 -coupling coefficient g' b3 /2tt — 46 Hz. In both the 
spring models of (a) and (b) the intra-cavity photon number 
versus detuning ti(Al) is calibrated from the known input 
laser power, cavity linewidth, and on-resonance transmission 
contrast. 


supermode resonance, which away from the anti-crossing 
point in this case is the slot-mode a\. In the sideband 
unresolved regime (u; m ft), the dynamic optical spring 
effect has a dispersive lineshape centered around the 
optical resonance frequency, with optical softening of the 
mechanical resonance occurring for red detuning (A l < 
0 ) and optical stiffening occurring for blue detuning 
(Al > 0 ). 

A fit to the measured frequency shift versus A^ is 
performed using the linear optomechanical coupling rate 
g ai ,b 3 as a fit parameter. The resulting optomechanical 
coupling rate which best fits the data is shown in 
Figure [ 6^1 as a red curve, and corresponds to g ai ,6 3 /27T — 
1.72 MHz. Using x zp f = 16 fm for the 63 mechanical 
mode, this corresponds to g ai ,b 3 /2n = 107 GHz/nm. 
Note that this is slightly smaller than the value measured 
indirectly from the dispersion in the anti-crossing curve of 
Fig- m however, that value relied on the simulated value 
for <a cap which is quite sensitive to the actual fabricated 
dimensions and stiffness of the structure. For the smallest 
splitting measured in this work (2J/2tt = 12 GHz), we 
get an estimated value for the x 2 -coupling to the 63 
mode from the dynamic optical spring measurements of 
g'b 3 / 27T = 245 Hz ( g'bJ’l-K = 0.96 THz/nm 2 ). 

An entirely different dynamics occurs at the 
anti-crossing point where x 2 optomechanical coupling 
dominates. Optical pumping of the supermode 


resonances near the anti-crossing point gives rise to an 
optical spring shift which depends on the static (i.e., not 
how it modulates with motion) value of the intracavity 
photon number. Due to the opposite sign of the quadratic 
dispersion of the upper and lower optical supermode 
branches, optical pumping of the upper branch resonance 
leads to a stiffening of the mechanical structure whereas 
optical pumping of the lower branch leads to a softening 
of the structure (231IH]. The measured frequency shift 
of the 63 mechanical resonance for optical pumping of 
the upper branch cavity supermode (the even a + mode 
in this case) is shown in Fig. [ 6 ]d for a voltage setting 
on the capacitor electrodes of Vi = 10.6 V and V 2 = 
1 V. This position is slightly below the exact center 
of the anti-crossing point of Fig. S* so as to allow 
weak linear transduction of the 63 resonance. A rather 
large supermode splitting of 2J/27T = 50 GHz was also 
chosen to ensure that only the even a + supermode is 
excited, and that the contribution to the optical trapping 
(anti-trapping) by the lower branch a_ resonance is 
negligible. 

As per eqs. © and the mechanical frequency 
shift is approximately given by Acu m (Az y ) « 2g' b3 n+(Ajr J ), 
where n+(A^) is the average intra-cavity photon number 
in the a + supermode. Fitting this model to the data 
measured in Fig. [fjj} yields a value of g' h J 2tt = 46 Hz. 
This is slightly lower than the 60 Hz value expected for 
a splitting of 2 J/ 27 T = 50 GHz and the linear coupling 
rate of g ai ,b 3 /2 tt = 1.72 MHz determined from the 
dynamical optical spring effect, but consistent with our 
slight detuning of the structure from the exact center of 
the anti-crossing. 


V. DISCUSSION 

The quasi two-dimensional photonic crystal 
architecture as presented here provides a means of 
realizing extremely large dispersive x 2 -coupling between 
light and mechanics. This is due to the ability to 
co-localize optical and acoustic waves in a common 
wavelength scale volume, resulting in inherently large 
linear optomechanical coupling. Combined with an 
ability to engineer the optical mode dispersion to 
allow for a tunable degree of optical mode splitting, 
the x 2 -coupling can be even further enhanced. It is 
interesting to consider then, just how far this technology 
could be pushed given recent technical advances made 
in the area of photonic crystals and optomechanical 
crystals. 

We consider here the feasibility of a QND measurement 
of phonon number, although similar parameters would 
enable a measurement of phonon shot noise [15], a QND 
measurement of photon number Et and the cooling 
and squeezing of mechanical motion [T 6 HT 8 ] . In the 
quasi-static limit as realized in this work (| J\ o; m ), the 
optical resonance shift per phonon is Acj = 2 g' = g 2 / J . 
If the lower frequency optical resonance (a_ in the case 
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J > 0 ) is used to probe the system, then roughly 
the photons emitted per unit time from the a_ cavity 
mode would change by (n_ft_)(Acj/ft_) upon a single 
phonon jump in the mechanical resonator. Assuming 
shot-noise limited detection over a measurement time r, 
the signal to noise ratio (SNR) for a phonon jump is given 
approximately by, 


(n_ft_) 2 (Ac<;/ft_) 2 T 2 
ri-K-T 


SNR: 


i-Auj 2 


ft_ 


( 12 ) 


The corresponding phonon jump measurement rate 
follows from the term in the bracket of eq. (12), r meas = 
[(A W )y._]n_ = [4(^) 2 /K_]n_. ^ 

This measurement rate should be compared against 
the decoherence rate of the mechanical resonator. The 
thermal decoherence rate is T t h = (n t h + l) 7 i, where 
nth is the Bose occupation factor depending on the bath 
temperature (T&), and 7 \ is the intrinsic mechanical 
damping rate to the bath. At X 5 = 4 K similar 
silicon photonic crystal devices have been operated with 
intra-cavity photon numbers of 10 3 and mechanical 
Q-factor as large as 7 x 10 5 DU- For the device studied 
here (cj m /27r ~ 10 MHz, g' /27r = 240 Hz, ft_/27r = 
5 GHz), the phonon jump measurement rate would be 
T m eas/27r ~ 46 mHz, while the thermal decoherence rate 
at T\) = 4 K and for Q m = 7 x 10 5 is T t h/27r ~ 125 kHz. 
Significant improvements in the measurement rate can be 
realized with improved optical Q-factor. Recent work by 
Sekoguchi, et al. [42], has shown that optical Q-factors 
of order 10 7 can be realized in similar planar 2D silicon 
photonic crystals in the telecom band, corresponding 
to a minimum cavity decay rate of k/2h = 20 MHz. 
By proper tuning of the double-slotted photonic crystal 
structure, the optical mode splitting 2 J could be reduced 
down to a minimum resolvable value equal to ft, yielding 
an x 2 -coupling value of g' /27r « 100 kHz and a phonon 
jump measurement rate of T meas /27r ~ 2 MHz. 

In order to realize a sideband-resolved system, higher 
mechanical resonant frequencies must also be employed. 
Numerical simulations [32 indicate that higher-order 
modes of the central nanobeam can maintain significant 
optomechanical coupling, with g/2 tt ~ 0.4 MHz for 
the seventh-order in-plane mechanical resonance at 
cj m /27r = 225 MHz. Tuning the structure such that the 
mode splitting is nearly resonant with the mechanical 
frequency, g <C \S = | 2 J| — u; m | cj m ,| 2 Jj, greatly 
enhances the frequency shift per phonon as per eq. 

Acj = g 2 /S. For similar cavity conditions as above 
(n_ = 10 3 , ft_/2?r = 20 MHz), and assuming S = 
10g, yields a measurement rate of r meas « 80 kHz. 
This is comparable to the thermal decoherence rate at 
Tb = 4 K assuming a similar mechanical Q-factor for 
these higher frequency modes. Recent measurements 
at bath temperatures of TJ, < 100 mK, however, have 
shown that mechanical Q-factors in excess of 10 7 can 
be realized in silicon using phononic bandgap acoustic 
shielding patterned in the perimeter of the device m- 


At these temperatures we can expect a bath occupancy 
of n t h ~ 10 , and with an acoustic bandgap shield, a much 
smaller thermal decoherence rate of T t h/27r ~ 300 Hz. A 
comparable measurement rate could then be employed 
with a much weaker optical probe corresponding to an 
intra-cavity photon number of n_ ~ 10 . 

The most challenging aspect of a QND phonon 
number measurement, however, is the optically induced 
mechanical decay due to residual back-action stemming 
from the linear (in x) cross-coupling of the cavity 
supermodes mini. This parasitic back-action damping 
of the mechanical resonator occurs through a process, 
for example, in which a photon is scattered from the 
driven a_ mode into the a + mode where it decays 
into the optical bath, absorbing a phonon in the 
process. The optically-induced mechanical decaj rate for 
a n^-phonon Fock state is given by r opt ~ (g/S) n&n_ft_ 
(= ((^') 2 /|2J|) n^n_ft + in the quasi-static limit) [44j. 
Comparing to the phonon jump measurement rate we 
see that only in the vacuum strong coupling limit (<?/ft > 
1) can one realize a continuous QND measurement of 
phonon number, 


^meas ^ P < f (13) 

r op t 77/5ft_|_ft_ \ft/ 

Note that a more careful analysis mm indicates that a 
limit of g > ft^ need only be met, where ft$ is the intrinsic 
damping of the optical cavity excluding loading of the 
cavity by measurement channels. A ratio of p/ft* ~ 0.007 
has previously been realized in silicon optomechanical 
crystals [44] . In the case of the double-slotted photonic 
crystal structure studied here, fabrication of nanoscale 
slots as small as s = 25 nm [45] would increase the linear 
optomechanical coupling between a± cavity supermodes 
to p/27r ^ 10 MHz. With this advance, and in 
conjunction with an increase of the optical Q-factor to 
10 7 [42], it does seem feasible in the near future to reach 
the vacuum strong coupling limit which would enable 
QND phononic and photonic measurements as proposed 
in Ref. [14] . 


ACKNOWLEDGMENTS 

The authors would like to thank M. Davanco and 
Max Ludwig for fruitful discussions. This work was 
supported by the AFOSR Hybrid Nanophotonics MURI, 
the Institute for Quantum Information and Matter, 
an NSF Physics Frontiers Center with support of the 
Gordon and Betty Moore Foundation, the Alexander 
von Humboldt Foundation, the Max Planck Society, 
and the Kavli Nanoscience Institute at Caltech. FM 
acknowledges support from the DARPA ORCHID 
program, ERC OPTOMECH, and ITN cQOM. TKP 
gratefully acknowledges support from the Swiss National 
Science Foundation. MK was supported by a College 









11 


Doctoral Fellowship from the University of Southern California. 


[1] M. Aspelmeyer, T. J. Kippenberg, and F. Marquardt, 
Rev. Mod. Phys. 86 , 1391 (2014) 

[2] T. P. Purdy, D. W. C. Brooks, T. Botter, N. Brahms, 
Z.-Y. Ma, and D. M. Stamper-Kurn, Phys. Rev. Lett. 
105, 133602 (2010) 

[3] M. Eichenfield, J. Chan, R. Camacho, K. Vahala, and 
O. Painter, Nature 462, 78 (2009) 

[4] B. Abbott et al. (LIGO Scientific Collaboration), New J. 
Phys. 11, 073032 (2009) 

[5] A. H. Safavi-Naeini, T. P. Mayer Alegre, J. Chan, 

M. Eichenfield, M. Winger, Q. Lin, J. T. Hill, D. E. 
Chang, and O. Painter, Nature 472, 69 (2011) 

[ 6 ] S. Weis, R. Riviere, S. Deleglise, E. Gavartin, O. Arcizet, 
A. Schliesser, and T. J. Kippenberg, Science 330, 1520 
( 2010 ), 

[7] A. D. O’Connell, M. Hofheinz, M. Ansmann, R. C. 
Bialczak, M. Lenander, E. Lucero, M. Neeley, D. Sank, 
H. Wang, M. Weides, J. Wenner, J. M. Martinis, and 
A. N. Cleland, Nature 464, 697 (2010) 

[ 8 ] J. Chan, T. P. M. Alegre, A. H. Safavi-Naeini, J. T. 
Hill, A. Krause, S. Groblacher, M. Aspelmeyer, and 
O. Painter, Nature 478, 89 (2011) 

[9] J. D. Teufel, T. Donner, D. Li, J. W. Harlow, M. S. 
Allman, K. Cicak, A. J. Sirois, J. D. Whittaker, K. W. 
Lehnert, and R. W. Simmonds, Nature 475, 359 (2011), 

[10] D. W. C. Brooks, T. Botter, S. Schreppler, T. P. Purdy, 

N. Brahms, and D. M. Stamper-Kurn, Nature 488, 476 

( 2012 ), 

[11] A. H. Safavi-Naeini, S. Groblacher, J. T. Hill, J. Chan, 
M. Aspelmeyer, and O. Painter, Nature 500, 185 (2013), 

[12] J. D. Thompson, B. M. Zwickl, A. M. Jayich, 
F. Marquardt, S. M. Girvin, and J. G. E. Harris, Nature 
452, 72-75 (2008) 

[13] H. Miao, S. Danilishin, T. Corbitt, and Y. Chen, Phys. 
Rev. Lett. 103, 100402 (2009) 

[14] M. Ludwig, A. H. Safavi-Naeini, O. Painter, and 
F. Marquardt, Phys. Rev. Lett. 109, 063601 (2012) 

[15] A. A. Clerk, F. Marquardt, and J. G. E. Harris, Phys. 
Rev. Lett. 104, 213603 (2010) 

[16] M. Bhattacharya, H. Uys, and P. Meystre, Phys. Rev. 
A 77, 033819 (2008) 

[17] A. Nunnenkamp, K. Bqrkje, J. G. E. Harris, and S. M. 
Girvin, Phys. Rev. A 82, 021806 (2010) 

[18] M. Asjad, G. S. Agarwal, M. S. Kim, P. Tombesi, G. D. 
Giuseppe, and D. Vitali, Phys. Rev. A 89, 023849 (2014), 

[19] M. R. Vanner, Phys. Rev. X 1, 021011 (2011). 

[20] A. H. Safavi-Naeini, Quantum Optomechanics with 
Silicon Nanostructures , Ph.D. thesis, California Institute 
of Technology (2013). 

[21] H. Kaviani, C. Healey, M. Wu, R. Ghobadi, A. Hryciw, 
and P. E. Barclay, Optica 2 , 271 (2014). 

[22] J. T. Hill, Nonlinear optics and wavelength translation via 
cavity-optomechanics, Ph.D. thesis, California Institute 
of Technology (2013). 

[23] C. Doolin, B. D. Hauer, P. H. Kim, A. J. R. MacDonald, 


H. Ramp, and J. P. Davis, Phys. Rev. A 89, 053838 
(2014), 

[24] G. A. Brawley, M. R. Vanner, P. E. Larsen, S. Schmid, 
A. Boisen, and W. P. Bowen, arXiv: 1404.5746 (2014) 

[25] N. E. Flowers-Jacobs, S. W. Hoch, J. C. Sankey, 
A. Kashkanova, A. M. Jayich, C. Deutsch, J. Reichel, 
and J. G. E. Harris, Appl. Phys. Lett. 101, 221109 
( 2012 ) 

[26] J. C. Sankey, C. Yang, B. M. Zwickl, A. M. Jayich, and 
J. G. E. Harris, Nature Phys. 6 , 707-712 (2010) 

[27] D. H. Santamore, A. C. Doherty, and M. C. Cross, Phys. 
Rev. B 70, 144301 (2004). 

[28] A. M. Jayich, J. C. Sankey, B. M. Zwickl, C. Yang, J. D. 
Thompson, S. M. Girvin, A. A. Clerk, F. Marquardt, 
and J. G. E. Harris, New J. Phys. 10, 095008 (2008). 

[29] A. A. Gangat, T. M. Stace, and G. J. Milburn, New J. 
Phys. 13, 043024 (2011). 

[30] K. Stannigel, P. Komar, S. J. M. Habraken, S. D. 
Bennett, M. D. Lukin, P. Zoller, and P. Rabl, Phys. 
Rev. Lett. 109, 013603 (2012) 

[31] M. Ludwig, Collective quantum effects in optomechanical 
systems , Ph.D. thesis, University of Erlangen-Nuremberg 
(2013). 

[32] M. Kalaee, T. K. Parai’so, and O. Painter, “Design of 
a quasi- 2 d photonic crystal optomechanical cavity with 
tunable, large x 2 -coupling,” (2015). 

[33] S. G. Johnson and J. D. Joannopoulos, Opt. Express 8 , 
173 (2001) 

[34] The MIT Photonic-Bands (MPB) software package, 
http://ab-initio.mit.edu/mpb. 

[35] COMSOL Multiphysics 3.5, http://www.comsol.com/. 

[36] A. H. Safavi-Naeini, T. P. M. Alegre, M. Winger, and 

O. Painter, Appl. Phys. Lett. 97, 181106 (2010) 

[37] M. Winger, T. D. Blasius, T. P. M. Alegre, A. H. 
Safavi-Naeini, S. Meenehan, J. Cohen, S. Stobbe, and 
O. Painter, Opt. Express 19, 24905 (2011) 

[38] C. P. Michael, M. Borselli, T. J. Johnson, C. Chrystal, 
and O. Painter, Opt. Express 15, 4745 (2007) 

[39] M. Eichenfield, R. Camacho, J. Chan, K. J. Vahala, and 
O. Painter, Nature 459, 550 (2009) 

[401 M. Bao and H. Yang, Sensors and Actuators A-Physical 
136, 3 (2007). 

[41] D. Lee, M. Underwood, D. Mason, A. B. Shkarin, S. W. 
Hoch, and J. G. E. Harris, Nat. Commun. 6 , 6232 (2015) 

[42] H. Sekoguchi, Y. Takahashi, T. Asano, and S. Noda, 
Opt. Express 22, 916 (2014) 

[43] S. M. Meenehan, J. D. Cohen, G. S. MacCabe, F. Marsili, 
M. D. Shaw, and O. Painter, arXiv:1503.05135 (2015). 

[44] J. Chan, A. H. Safavi-Naeini, J. T. Hill, S. Meenehan, 
and O. Painter, Appl. Phys. Lett. 101, 081115 (2012). 

[45] A. Pitanti, J. M. Fink, A. H. Safavi-Naeini, J. T. Hill, 
C. U. Lei, A. Tredicucci, and O. Painter, Opt. Express 
23, 3196 (2015). 














